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Abstract 

The fairly recent discovery of "quasicrystals" , whose X-ray diffrac- 
tion patterns reveal certain peculiar features which do not conform 
with spatial periodicity, has motivated studies of the wave-dynamical 
implications of "aperiodic order" . Within the context of the radiation 
properties of antenna arrays, an instructive novel (canonical) exam- 
ple of wave interactions with quasiperiodic order is illustrated here for 
one-dimensional (1-D) array configurations based on the "modified- 
Fibonacci" sequence, with utilization of a two-scale generalization of 
the standard Poisson summation formula for periodic arrays. This 
allows for a "quasi-Floquet" analytic parameterization of the radiated 
field, which provides instructive insights into some of the basic wave 
mechanisms associated with quasiperiodic order, highlighting similar- 
ities and differences with the periodic case. Examples are shown for 
quasiperiodic infinite and spatially-truncated arrays, with brief dis- 
cussion of computational issues and potential applications. 
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1 Introduction 

The recent discovery (1984) of "quasicrystals" , i.e., certain metallic alloys 
whose X-ray diffraction patterns contain bright spots displaying symme- 
tries (e.g., 5- fold) which are incompatible with spatial periodicity [T], [2], 
has stimulated a growing interest in the study of aperiodic order and its 
wave-dynamical properties. 

In electromagnetics (EM) engineering, use of random or deterministic 
aperiodic geometries has been customary within the framework of antenna 
array thinning [H] [S], whereas multiperiod configurations have recently been 
proposed for optimizing the passband/stopband characteristics of frequency 
selective surfaces jU] and photonic bandgap (PBG) devices j7]. In jS], we 
explored the radiation properties of two-dimensional (2-D) antenna arrays 
based on the concept of "aperiodic tiling" 0, [TU], which had previously 
found interesting applications in the field of PBG devices (see the brief sum- 
mary and references in [5|). 



In this paper, we turn our attention to the radiation properties of a dif- 
ferent category of aperiodic structures consisting of two-scale Fibonacci-type 
sequences Materials exhibiting this type of aperiodic order, technically 
called "quasiperiodicity" (see the definition in Sec. 12. 2|) . were first fabricated 
in 1985 as GaAs-AlAs heterostructures (1-D multilayers) [T2]; their techni- 
cal classification as "quasicrystals" is still debated ^3]. The wave-dynamical 
properties of 1-D and 2-D Fibonacci-type structures have been widely investi- 
gated, theoretically and experimentally, in quantum mechanics, acoustics and 
EM (see [H|-[2B] for a sparse sampling). Particularly interesting outcomes 
concern the self-similar fractal (Cantor-type) nature of the eigenspectra |14j . 
[To] . |2Uj . and the possible presence of bandgaps [T^j, |2Ij, omnidirectional 
reflection properties (221, [23] . and localization phenomena [Hj, [THj, [23] • 

Here, we concentrate on the study of the radiation properties of a simple 
class of 1-D antenna arrays based on the so-called "modified-Fibonacci" se- 
quence. This novel prototype array configuration appears particularly well 
suited to exploration of some basic characteristics of wave interactions with 
quasiperiodic order. First, this array model is amenable to analytic parame- 
terization via a generalized Poisson summation formula, by exploiting some 
recent results in (2ZI (this paper is easily accessible through the web). Note 
that, in principle, one can derive generalized Poisson summation formulas 
that accommodate a variety of "nonperiodicity scales" (see, e.g., |2HJ). Here, 
we consider only two scales, d\ and c?2- This opens up the possibility of 
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extending the Floquet-based parameterization of infinite and semi-infinite 
time-harmonic periodic arrays in |2H], [30j to the case of two-scale quasiperi- 
odic arrays. Next, the inherent degree of freedom in the choice of the ratio 
between the two scales can be used to study the "transition" from periodic 
(d\ = d 2 ) to quasiperiodic {d\ ^ d 2 ) order, so as to better understand the 
quasiperiodicity-induced footprints in the wave dynamics. Finally, although 
the main focus of this preliminary investigation is on wave-dynamical phe- 
nomenologies, computational and applicational issues for a test example are 
briefly addressed as well, including possible exploitation of the two-scale de- 
gree of freedom for pattern control in practical applications. 

The paper is organized as follows. In Section |2J the problem geometry is 
described, and the modified-Fibonacci sequence is introduced together with 
general aspects of quasiperiodicity. In Section El the generalized Poisson 
summation formula for modified-Fibonacci arrays is introduced, and its sim- 
ilarities and differences with the periodic case are discussed. In Section HI a 
"quasi-Floquet" (QF) parameterization of the radiated field for infinite and 
semi-infinite modified-Fibonacci arrays is derived, paralleling [201 • Nu- 
merical results and potential applications are illustrated in Section^ followed 
by brief concluding remarks in Section |3 

2 Background and Problem Formulation 

2.1 Geometry 

Referring to the geometry depicted in Fig. we begin with an infinite phased 
line array of z-directed electric dipoles, subject to uniform unit-amplitude 
time-harmonic exp(juot) excitation, described by the current distribution 

oo 

f( z )= $( z - z rn)exp(-jr]k z m ) , (1) 

m=— oo 

where k = Uy/e fi = 2tt/\ is the free-space wavenumber (with A being 
the wavelength), and — 1 < rj < 1 describes the inter-element phasing. The 
dipole sequence {z m }m=-oo' which is restricted to two possible inter-element 
spacings d\ and d 2 < d\ (see Fig. |TJ), is chosen according to the modified- 
Fibonacci rule, whose properties are summarized below. Spatial truncation 
effects will be discussed in Sec. 14.21 
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2.2 Modified-Fibonacci Sequence: An Example of 
Quasiperiodic Order 

The Fibonacci sequence jTTj, introduced in 1202 by the Italian mathemati- 
cian Leonardo da Pisa (Fibonacci) (ca. 1170-ca. 1240), in connection with 
a model for rabbit breeding, is probably the earliest and most thoroughly 
investigated deterministic aperiodic sequence. Since then, it has found ap- 
plications in many different fields, owing to its intimate relation with one of 
the most pervasive mathematical entities, the Golden Mean [11]. In its sim- 
plest version, the Fibonacci sequence can be generated from a two-symbol 
alphabet A = {a,b}, by iteratively applying the substitution rules 

a — > ab, b — > a, (2) 

so as to construct a sequence of symbolic strings s n 

b =>• a =>- ab =>- aba =>■ abaab =>■ abaababa =>■ abaababaabaab ■ ■ ■ (3) 

Note that the string at each iteration is obtained as the concatenation of 
the two preceding ones (s n = s„_i U s n _2). The process can be iterated ad 
infinitum, yielding an infinite sequence of "a" and "6" symbols which seems 
to display no apparent regularity, but actually hides a wealth of interesting 
properties (see jTTj for details). For instance, in the limit of an infinite 
sequence, it can be shown that the ratio between the numbers of "a" and 
"6" symbols (N a and iV&, respectively) approaches the Golden Mean r = 

(i+V5)/2im, 

hm§ = r = l±^, N a + N b = N. (4) 

In connection with the antenna array problem of interest here, there are 
several ways of embedding the above-introduced Fibonacci-type aperiodic 
sequence. One possibility would be keeping the geometry periodic (i.e., uni- 
form inter-element spacing) and associating with the "a" and "6" symbols in 
the Fibonacci sequence two possible current amplitudes. Another possibil- 
ity, pursued in this investigation, assumes uniform excitation and associates 
the "a" and "6" symbols with two possible inter-element spacings, d\ and 
d-2 < d\, respectively (see Fig. [TJ. Accordingly, the dipole positions z m in 
(JTJ) can be obtained from the symbolic strings in (J3J) or, in an equivalent and 
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more direct fashion, via |2*7j 

where ||-|| denotes the nearest-integer function 

m, m < x < m + \ 
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m + 1, m + | < x < m + 1. ^ 

Following the particular case <i 2 = rfi/V in (0) will be referred to as 
"standard Fibonacci". The general case, in which the scale ratio v = d 2 jd\ 
is left as a degree of freedom, will be referred to as "modified Fibonacci", 
and represents one of the simplest extensions of the Fibonacci sequence. The 
reader is referred to [2H], PI] for other examples of possible exten- 

sions/generalizations. 

Note that the modified-Fibonacci sequence in (JHJ) includes as a limit the 
•periodic case [d\ = d 2 ). It can be proved that, with the exception of this de- 
generate case, the sequence in (jSJ) and the corresponding dipole distribution 
in Q display a "quasiperiodic" character, which represents one of the most 
common and best known forms of aperiodic order. The concept of quasiperi- 
odicity stems from the theory of "almost-periodic" functions developed by 
H. Bohr during the 1920s Ell- in essence, a quasiperiodic function 
can be uniformly approximated by a generalized Fourier series containing a 
countable infinity of pairwise incommensurate frequencies generated from a 
finite-dimensional basis (see [HH], |M] f° r more details). 

It follows from Q that, in the infinite-sequence limit, the average inter- 
element spacing d av is j^Zj 

rdi + d 2 

d av = ~ ~ ■ (7) 
1 + T 

Anticipating the analytic derivations and parametric analysis in Sees. Eland 
it is expedient to parameterize the sequence in © in terms of the average 
spacing d av in ((7|) and the scale ratio v = d 2 /di, rewriting d\ and d 2 as 

di = [ 1 + T \ d av , d 2 = udi, < v < 1. (8) 
\v + V 

The modified-Fibonacci sequence in © admits an instructive alternative 
interpretation in terms of a "cut-and-project" graphic construction |2Zj, as il- 
lustrated in Fig. |21 Cut-and-project schemes are systematic tools for generat- 
ing quasiperiodic sets via projection from higher- dimensional periodic lattices 
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(see [TU] for details). In our case, one starts from a 2-D square lattice of side 
do = divl + v 2 traversed by a straight line with slope tan0o = 1/t. Those 
lattice points whose "Voronoi cell" [TU] (light-shaded square cell of side do 
centered around the point, in Fig. |2J) is crossed by the line (or, equivalently, 
those falling within the rectangular window of size h = d (l + t)/\/2 +t 
centered around the line, in Fig. are orthogonally projected onto another 
straight line (z-axis in Fig. |2J) with slope tan0 = v to yield the desired 
modified-Fibonacci sequence in (JHJ. For the standard-Fibonacci sequence 
[y = 1/t), the two lines coincide (0 = <po) and the above scheme becomes 
equivalent to the canonical cut-and-project scheme described in jTUj. 



3 Generalized Poisson Summation Formula 

For periodic structures, Floquet theory provides a rigorous and powerful 
framework for spectral- or spatial-domain analytic and numerical analysis. 
In this connection, the Poisson summation formula |3S] can be utilized to sys- 
tematically recast field observables as superpositions of either individual or 
collective contributions. Problem-matched extensions have also been devel- 
oped to accommodate typical departures from perfect periodicity in realistic 
structures, such as truncation (finiteness) and smooth perturbations in the 
spatial period as well as in the excitation (tapering) Po]-[3$]. 

Considering the degenerate periodic limiting case (di = d 2 = d, i.e., z m = 
md) of the modified-Fibonacci sequence in (J3J, the corresponding periodic 
limit /( per ) of the current distribution in (Q) can be recast, via the standard 
Poisson summation formula [33], as 

oo 

f( per \z) = >J S(z — md) exp (—jrjk md) (9a) 

m=— oo 

= - ^ exp (-jk zg z) , k zq = k r] + — . (9b) 

q= — oo 

The m-indexed individual dipole contributions in ()9a|) are thereby recast into 
the infinite superposition of linearly smoothly-phased g-indexed equivalent 
line source distributions in (|9b|). Remarkably, a similar re-parameterization 
is always possible for the general (i.e., quasiperiodic) case of the modified- 
Fibonacci sequence in (J3J), since it can be shown that the spatial Fourier 
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transform (plane-wave spectrum) of (£[]) can be written as j2Zj (see also the 
Appendix) 

/oo 
f(z)exp(jk z z)dz 
-oo 
oo 

= 22 e wUk z z m )exp(-jr]k z m ) (10a) 

m=— oo 

= -3 S qiq2 5(k z — k zqiq2 ), (10b) 



91,92=— oo 



where d au is defined in flj|), and the amplitude coefficients S qiq2 and the 

9192 

sin W n 



'qiqi 

spatial frequencies k zqiq2 are given by 



_ ■ ■ qiQ'2 

9192 ~ W 

vv qiq2 



\r 71 1 a a\ 7r ( 1 + r )(9i -92^) ..... , 

^9192 = (<?l rf l - ?2^ 2 ) = , (11) 

d av V + T 



kzq!q 2 koTj -\- (X qiq2 , Ot qiq2 j f_ , t\ ' (1^) 



2ir (q 1 + g 2 r 
d a „ (r + 1) 

It then follows via straightforward inverse Fourier transform of (|10b|) that 



1 x 

f( z ) = l~ S ng2^p(-jk Z q iq2 z) . (13) 



dav 

<?i,<?2=-oo 



Equation (JT3J) thus represents a generalization of the standard Poisson sum- 
mation formula for periodic arrays in Q to the more general quasiperiodic 
modified-Fibonacci array in ((HJ), and will accordingly be referred to as gen- 
eralized Poisson summation formula. It is readily verified that for the special 
case of periodic arrays {d\ — d% = d, i.e., v = 1), one obtains 

d av = d, W qiq2 = (qi - q 2 )ir, S qiq2 = 5 qiq2 , (14) 

with 5 qiq2 denoting the Kronecker delta, and thus (fTSjl reduces to the standard 
Poisson summation formula in (j9b|) . 

A few observations are in order to highlight similarities and differences 
between and (JOJ). Similar to (JSJ), the individual dipole contributions in 
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(JHJ) are recast via (JTHj) as a superposition of collective contributions arising 
from linearly smoothly-phased equivalent line source distributions. However, 
at variance with the (/-indexed single infinity of equispaced spatial frequen- 
cies in (j9b|) . the generalized Poisson summation formula in ()13j) entails a 
(gi, ^-indexed double infinity of generally pairwise-incommensurate spatial 
frequencies. This could be expected recalling the definition of the quasiperi- 
odic functions given in Sec. 12.21 Note that the spatial frequencies k zqiq2 in 
(H2J) depend on the average inter-element spacing d av , but are independent 
of the scale ratio v = g^/^i- The dependence on v shows up (via W qiq2 ) in 
the amplitude coefficients S qiq2 in (fTTj) . which are not constant as for the pe- 
riodic case in (|9b|). It is worth pointing out that such dependence is smooth, 
and consequently no abrupt transition occurs when v is varied between ratio- 
nal and irrational values, i.e., going from commensurate to incommensurate 
scales d\ and d 2 - However, for case of commensurate scales, 

" = T = -i Pi.^N, ( 15 ) 
one can readily verify that 

kz(q 1 +Mp 1 )(q 2 +Mp 2 ) 
W( qi+Mpi )(q 2 +Mp2) 

and thus the array spectrum in (jlObjl is periodic with period 2717)2/^1 • It can 
be verified that the opposite implication is also true, and thus the spectrum 
is periodic if and only if the scales di and d 2 are commensurate, although 
the array remains aperiodic in physical space |27j . 

4 Radiated Field 

In order to further explore possible similarities and differences between the 
wave phenomenologies associated with periodicity and quasiperiodicity, it is 
instructive to consider the radiated fields. To this end, in what follows, we 
extend the Floquet-based frequency-domain analysis in [2H], [30], for infi- 
nite and semi-infinite periodic phased arrays of dipoles, to the quasiperiodic 
modified-Fibonacci case. As in [2H|, [HO], attention is restricted to the z- 
directed vector potential A(r) = A(r)u z , with r = (z, p) and with u 2 denot- 
ing a z-directed unit vector, from which all field quantities of interest can be 
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computed. Whenever applicable, partial results from [2H|, [20] are recalled 
and used, without going into the details of the technical derivations. 

4.1 Infinite Arrays 

For infinite arrays, proceding analogous to [21] , the generalized Poisson sum- 
mation formula in (J13|) can be used to recast the element-by-element (spher- 
ical wave) synthesis of the field potential, 

At \ f°° ft >\ eic P{-]hR), , ^ . , , 

A(r) = J f(z) dz = Y: A m (r) 

J ~°° m=-oo 

= Y 6XP 4 ^° Rm exp (-jk r}z m ) , (17) 

m=-oo m 

where R = p 2 + (z — z') 2 and R m = a/ p 2 + (z — z m ) 2 , in terms of a "quasi- 
Floquet" (QF) representation 

oo 

A(r)= ( 18 ) 

91)92=— oo 



< 2 (r) 



exp (-jkoR) 



dn 



exp (-jk 



AtxR 



exp (-jk zqig2 z') dz' 



Ajd, 



Hq ^ (kpq lQ2 p) . 



(19) 



In (119)1 . Hq 2 ' ) denotes the zeroth-order Hankel function of the second kind 

fcL„„, lm(k pqiq2 ) < de- 



(line-source Green's function), and k pqig2 = y k% 

note the QF wave radial wavenumbers related to the ^-domain wavenumbers 
k zqiq2 in (fT2~|l. Recalling the asymptotic (|fc pgig2 p| — > oo) expansion of the 
Hankel function, one obtains 

< 2 « 



d av ^Kk 



P9192/ 



x exp 



J ( k pqiq2 p + 



7T 



(20) 
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from which it is recognized that 



k 



, ^ Qi , fd av 
< k 0l q 2 ^ ± -T- 



) 



(lT^)(l + r) 



r 



(21) 



corresponds to radially-propagating QF waves, whereas |/c 2?ig2 | > k corre- 
sponds to radially-evanescent QF waves. Accordingly, sufficiently far from 
the array axis, the potential field in (|18j) will be synthesized in terms of 
conical propagating QF waves (heavy solid arrow in Fig. |2J), with arrival 
directions 



In the (gi, g 2 )-plane, the propagating spectral range in (|2Tj) is mapped into 
an infinite strip, thus indicating the general presence of an infinite number of 
propagating QF waves. This is illustrated in Fig. HJ for a nonphased (rj = 0) 
array with d av = l.lAo and various values of the scale ratio v = dz/dx. 
Also shown, for comparison, is the periodic case (Fig. EJa)), which entails 
a finite number of propagating Floquet waves. In the general quasiperiodic 
case, it is observed that the propagating spectral range is vastly populated, 
irrespective of the commensurate (Fig. EJb)) or incommensurate (Fig. EJc)) 
character of the scales d\ and d<i- However, it follows from (jllj) that moving 
toward large values of q\,q 2 (with opposite sign), the amplitude coefficients 
S qiq2 decay non-monotonically as ~ (gi — q 2 v)^ x ■ For a better quantitative 
understanding, Fig. |S] shows the direct mapping S qiQ2 vs. k zqiq2 for the 
same array configurations as in Fig. and \qi\, \q 2 \ < 50. Apart from the 
trivial periodic case (Fig. EJa)), one observes highly-populated spectra, which 
display perfect periodicity in the case of commensurate scales (Fig. E(b)) and 
only some loose repetitiveness otherwise (Fig. Efc)). In both vast 
majority of the spatial frequencies have amplitude coefficients significantly 
smaller (< — 20dB) than the dominant ones. 

4.2 Truncation Effects: Semi-Infinite Arrays 

Capitalizing on some analogies with the problem expounded in j3U], we have 
also studied truncation-induced diffraction effects. For the semi-infinite (m > 




(22) 
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0) version of the array in (Q), one obtains a truncated QF wave synthesis, 

oo 

A(r) = £> m (r) 



m=0 

oo 



Eexp (-jk R m ) 
— exp (-jkoVZr, 
ATTKm 



m=0 

oo 



Ai,{T> + £ (23) 



9i>92=-°o 

where ^ 192 are truncated QF wave propagators, 

lav JO 



A qiq2 (r) = exvHk zqiq2 z')dz> 

(2) 



exp ( jk zqiq . 2 z) f°° H (k pqiq2 p) 



8ird n 



r H ° ^ pqiq2p) dk z . (24) 



Exploiting the uniform high-frequency asymptotic approximation given in 
Eqs. (28)-(35) of [3U] for the spectral integral in (j2U), one obtains 

<*M - a&Wu 6C - &) + K«M ( 2 5) 

In fl23j) . is the QF wave propagator in (J2UJ), £/ denotes the Heavi- 

side unit-step function, P q ^ 2 delimits the shadow-boundary of each QF wave 
(for propagating QF waves, P q ^ 2 = /3 qig2 ). Moreover, A^ iq2 represents the 
diffracted QF wave emanating from the array tip, 

d ^ exp (-jk R d ) F^J 

91921 ; j47rd av k R d (cos (3 qiq2 - cos [3 d y 1 ; 

where R d = a/p 2 + z 2 , Pd = arctan(p/z), 7 9l92 = y/2k Q R d sin[(/3 9l92 - A*)/2], 
and F denotes the standard transition function of the uniform theory of 
diffraction (UTD), 



/"OO 

F(x) = 2jVxexp(jx) / exp (-j£ 2 ) d£ 



--7T < axg(x) < |. (27) 



11 



The wave phenomenologies are illustrated in Fig. |2l for the case of prop- 
agating (\k zqiq2 \ < ko) QF waves. The region of validity of each QF wave 
A®^ 2 (heavy arrow in Fig. |3J) arriving from direction f3 qiq2 is now limited 
by a conical shadow boundary. The spherical tip-diffracted QF wave A^ iq2 
(dashed arrow in Fig. |3J) arriving from direction j3d ensures, via the transition 
function in (J27j) . continuity of the wavefield across the parabolic transition 
region (gray shading in Fig. ^ surrounding the shadow boundary cone. 

Note that evanescent (\k zqiq2 \ > ko) QF waves yield negligible contribu- 
tions at observation points far from the array axis, but they excite propa- 
gating diffracted fields that need to be taken into account. As in |3U] (see 
the discussion after Eq. (35) there), these contributions are approximated 
via nonuniform asymptotics (F ~ 1 in (|26|) ) . We point out that, at variance 
with the periodic array case (cf. Eq. (36) in [30J ) , it is not possible here to 
recast the total spherical wave diffracted field in a more manageable form. 

5 Results and Potential Applications 

Although our main interest in this preliminary investigation is focused on 
wave-dynamical phenomenologies associated with radiation from quasiperi- 
odic antenna arrays, we briefly discuss some computational and applicational 
aspects with the hope of providing further insights. We stress that no at- 
tempt has been made at this stage to devise optimal computational schemes, 
nor do we deal with actual fabrication-oriented issues (feeding, matching, 
inter-element coupling, etc.). 

5.1 Numerical Results 

From a computational viewpoint, the actual utility of the QF syntheses 
in Sec. 0] (double summations involving an infinite number of propagating 
waves) might appear questionable as compared with brute-force element-by- 
element synthesis. However, as observed in Sec. 14.11 fsee also Fig. |HJ), under 
appropriate conditions, a large number of propagating QF waves could be 
weakly excited, thus suggesting the possibility, yet to be explored, of devising 
effective truncation schemes. Here, we try to quantify some of these aspects 
via numerical examples. We begin by considering a 101-element nonphased 
standard-Fibonacci array (rj = 0, u = 1/t, \m\ < 50 in (£Q)). The truncated 
QF synthesis developed in Sec. 14.21 for a semi-infinite array can readily be 
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exploited here by expressing the finite array interval as the difference be- 
tween two overlapping semi-infinite intervals. Figure El shows the near-field 
(R = lOOAo) QF synthesis results for d av = 0.5Ao, using the crudest possible 
truncation criterion based on retaining N p dominant propagating QF waves 
in (j2HJ). Also shown, as a reference solution, is the element-by-element syn- 
thesis in (II 7j) (with |m| < 50). It is observed that the QF synthesis with 
N p — 9 propagating waves provides reasonably good agreement, and even 
N p = 3 is still capable of fleshing out most of the wavefield structure. Sim- 
ilar statements can be made for the far-field pattern in Fig. [7| In order 
to better quantify the accuracy, and address convergence issues, we have 
computed the r.m.s. error 



where the subscripts and u Q Fn denote the reference solution and the QF 
synthesis, respectively. Figure |H1 shows the error behavior vs. the number N p 
of dominant propagating QF waves retained, for near-field synthesis (results 
for the far- field are practically identical). Also shown, for comparison, are 
results obtained by retaining a number N e = N p of dominant evanescent QF 
diffracted waves. It is observed that a moderate number (~ 10) of QF waves 
is capable of providing acceptable accuracy (A A ~ — 20dB). The contribu- 
tion of N e retained evanescent QF diffracted waves is practically negligible 
for N e < 20, but can give observable improvements at larger iV e -values. Intu- 
itively, one would expect the convergence behavior to improve in the presence 
of weaker aperiodicity {y « 1) and smaller average inter-element spacing, and 
viceversa. This is confirmed by the results shown in Fig. for two values 
of the inter-element spacing (d av = 0.5A and 0.75A ) and three values of the 
scale ratio [y = 0.9,0.75 and 1/r). 

To sum up, the above results seem to indicate that moderate-size QF 
syntheses, truncated by using even crude criteria, can still be capable of 
capturing the dominant features of the relevant wave dynamics. However, 
besides the computational convenience, which can become questionable if 
highly accurate results are needed, the QF parameterization can offer valu- 
able insights for judicious exploitation of the inherent degree of freedom (scale 
ratio v = d 2 /di) in the array (see, e.g., Sec. 15. 2|) . 




p_l) 2 \A^{R,e)-AQp{R,e)\ 2 de 



(28) 
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5.2 Potential Applications 

In view the discrete character of its spectrum (see (|10bjl ). the modified- 
Fibonacci array in (0) does not offer particular advantages within the frame- 
work of array thinning, as compared with periodic arrays. In this connection, 
other 1-D aperiodic sequences, such as the Rudin-Shapiro [20| (characterized 
by continuous spectra), might be worth being explored. 

However, the degree of freedom available in the choice of the scale ratio 
v = dijd\ can be exploited in principle to control the spectral properties and 
achieve, for instance, a multibeam radiation pattern. As a simple example, 
we consider a nonphased (rj = 0) configuration with average inter-element 
spacing d av < A . It is readily observed from (jTTjl . ()12|) that, besides the 
unit-amplitude QF wave (q x = q 2 = 0) (main beam at broadside), the only 
other propagating QF waves (secondary beams) are 

n(l + r)(l + u) 

qi = -q 2 = ±l, W qiq2 =- f ~>ir, 29a 

V + T 

gi = ±l,g 2 = 0, [W gig2 ]= ^ >ir, (29b) 

v + r 

gi = 0,g 2 = ±l, \W qiq2 \ = nU ^ T) <tt. (29c) 

v + r 



Note that, since |W gig2 | > n, the amplitude of the secondary beams in (j29aj) 
and (j29b|) will always be at least 13dB below the main beam, and thus at the 
same level of sidelobes as in finite-size periodic arrays, irrespective of the scale 
ratio v. Conversely, for the secondary beam in (J2~9"j) . one has |W gig2 | < n, and 
thus its amplitude can be controlled over a wide range by varying the scale 
ratio v. Restricting attention to the (q\ = 0, q 2 = 1) beam, one finds from 
(H2J) the corresponding direction /3 i (from endfire) and amplitude S i 

to be 



/3 i = arccos 



X t 



d av (l + t) 



[y + r) sin 
5*01 = — 



ttv(1+t) 

V+T 



7TZ/(1 + r) 



(30) 



The direction /3 i can thus be steered, by varying d av /\ , up to a maximum 
value (3 < ^ iax ' > = arccos [r/(l + r)] ~ 51.83°, corresponding to the maximum 
spacing (d av = A ) allowable to prevent emergence, in the visible range, of 
higher-order grating lobes. The amplitude Sqi can be controlled (from to 
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1, in principle) by varying the scale ratio v = d 2 jd\ (from 1 to 0). One thus 
obtains a multibeam radiation pattern with the possibility of controlling the 
secondary beam amplitude. Results for the (qi = 0,q 2 = —1) beam follow 
from symmetry considerations. Obvious fabrication-related issues prevent v 
from being exceedingly small, but values of v m 0.25 are sufficient to achieve 
secondary beam amplitudes ~ — 2dB. Figure ITU1 shows the radiation pattern 
(array factor) of a 101-element array, with d av /X and v chosen from (|3U|) so 
as to create a secondary beam at 45° with various amplitudes. It is observed 
that actual amplitude values are very close to the infinite-array predictions 
in (f50|). and that minor sidelobes never exceed the periodic-array sidelobe 
level (-13dB). 

It is hoped that the above observations might open up new perspec- 
tives for reconfigurable arrays. Note that rational values of v (commen- 
surate scales) correspond to configurations interpretable as periodic arrays 
with aperiodically- distributed lacunas. In these cases, one can think of easily- 
implementable reconfigurable strategies, based on switching from single-beam 
(periodic array) to multibeam (modified-Fibonacci) radiation patterns, via 
on-off selection of a set of antenna elements, keeping the average inter-element 
spacing d av < Ao- 

6 Conclusions and Perspectives 

A simple illustrative example of wave interaction with quasiperiodic order 
has been discussed in connection with the radiation properties of 1-D an- 
tenna arrays, utilizing the modified-Fibonacci sequence. A "quasi-Floquet" 
analytic parameterization, based on a generalized Poisson summation for- 
mula, has been derived for infinite and semi-infinite arrays. Computational 
issues and potential applications have been briefly addressed. 

It is hoped that the prototype study in this paper, through its instructive 
insights into some of the basic mechanisms governing wave interactions with 
quasiperiodic order, may lead to new applications in array radiation pattern 
control, in view of the additional degrees of freedom available in aperiodic 
structures. Accordingly, we are planning current and future investigations of 
modified-Fibonacci arrays that will emphasize the effects of the scale ratio 
parameter [y = d 2 /d\) on the input impedance as well as on the coupling to 
possible dielectric-substrate-induced leaky modes. Exploration of the radia- 
tion properties of antenna arrays based on other well-understood aperiodic 



15 



sequences (e.g., Thue- Morse, period-doubling, Rudin-Shapiro [20 j) is also 
being pursued. 



Appendix 

Pertaining to (TTO-(TT5h 



In [2I|, the modified-Fibonacci array spectrum in (fTU|) is computed using 
two equivalent approaches, one directly based on the cut-and-project scheme 
described in Sec. 12.21 (see also Fig. |2J), and the other based on an "average 
unit cell" method [32] • Although the two approaches are relatively simple in 
principle, their implementation is rather involved and not reported here for 
brevity. The reader is referred to j27], [32] for theoretical foundations and 
technical details, and to [10], [S] for an alternative approach applicable to 
rather general substitutional sequences. We point out that in [21] , the result 
for the radiation spectrum is given in normalized form, with a non-explicit 
multiplicative constant, and assuming zero phasing (rj = 0). In the paper 
here, the possible presence of phasing is accounted for via the spectral shift 
kof) in (jl2J) . The calculation of the proper multiplicative constant in ()10|) 
(and hence that in (fT3|) ) has been accomplished by first observing that for 
the modified-Fibonacci sequence in (0), the number of array elements falling 
within a window of width 2( approaches the ratio between the window width 
and the average inter-element spacing d av in (J7J), as ( — > 00, i.e., 

[U(z m + C)-U(z m -0]^-r, (31) 

m=—oo 

with U denoting the Heaviside step- function. It then follows from ([T]), as- 
suming 7] = and recalling ()31|). that 

1 [< 
lim — / f(z)dz 
C^oo ZQ J_t* 

1 °° 

= , lim ^7 E [U(z m + ()-U(z m -()} 

m=— 00 

(32) 



C^oo 2^ \d av J d a 

The proper multiplicative constants in (|T0|) and (|T3|l have accordingly 
been computed by requiring that (fT3|) satisfy (JH2]) for 77 = 0. 
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Figure 1: Problem schematic: An infinite (or semi-infinite) phased line array 
of ^-directed electric dipoles is considered. The dipole distribution z m , which 
features only two possible inter-element spacings d\ and d 2 < di, is chosen 
according to the modified-Fibonacci sequence in Also shown are the 
(z, p) and (R, 0) coordinate systems utilized. 



21 




Figure 2: "Cut-arid-project" graphic interpretation of the modified- Fibonacci 
sequence in ©. A 2-D square lattice of side d = d±\/l + v 2 is traversed 
by a straight line with slope tan0 o = 1/r, and those lattice points whose 
"Voronoi cell" (light-shaded square cell of side do centered around the point) 
is crossed by the line (or, equivalently, those falling within the dark-shaded 
rectangular window of size h = g? (1+t) / \/2 + r centered around the line) are 
orthogonally projected onto another straight line (z-axis) with slope tan0 = 
v. For the standard-Fibonacci sequence [y = 1/V), the two lines coincide 

(0 = 00). 
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Figure 3: Quasi-Floquet (QF) wave phenomenologies for the infinite and 
semi-infinite array (planar cut through conical wavefields). For the infinite 
array, the wavefield away from the array axis is synthesized in terms of prop- 
agating QF waves (heavy solid arrow) arriving from direction (3 qiq2 in 
(I22j) . For the semi- infinite array, the region of validity of each propagating 
QF wave is limited by a conical shadow boundary SBq 1(?2 with the same angle 
P qiq2 - A spherical tip-diffracted QF wave A qiq2 (dashed arrow) arriving from 
direction (3d ensures, via the transition function in (J27|) . continuity of the 
wavefield across the parabolic transition region (gray shading) surrounding 
the shadow boundary cone. 
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Figure 4: Modified-Fibonacci array with d a „ = 1.1A , 77 = 0, and various val- 
ues of the scale ratio v = d 2 /di. (qi, g 2 )-plane mapping of propagating (dark 
dots) and evanescent (light dots) amplitude coefficients S qiQ2 in (fTTj). Only 
\<li\, \<l2\ < 10 spatial frequencies are considered, with dot size proportional to 
coefficient amplitude. The two oblique dashed lines delimit the propagating 
spectral range in (|2Tj) . (a): u = 1 (periodic); (b): u = 0.75; (c): v — 1/r 
(standard-Fibonacci) . 
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Figure 5: As in Fig. H but \S qiQ2 \ vs. k zqiq2 , for |gi|, \q 2 \ < 50 and k zqiq2 < 7k . 
Due to symmetry, only positive spatial frequencies are shown. 
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e [deg] 

Figure 6: 101-element standard-Fibonacci {y — 1/r) array with d av = 0.5Ao 
and 7] = 0. Near-zone (R = lOOAo) normalized potential field scan. — 

— Reference solution (element-by-element summation); QF synthesis 

retaining N p = 3 dominant propagating waves (selected within \qi\, |g 2 | < 

50); QF synthesis retaining N p — 9 dominant propagating QF waves. 

Due to symmetry only positive angles are shown. 




6 [deg] 

Figure 7: As in Fig. U but far-field pattern. 
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Figure 8: Parameters as in Fig. |Hl but r.m.s. error AA in (|28|) vs. number N p 
of propagating waves retained in the near field (R = lOOAo) QF synthesis. 
Circular bullets: No evanescent diffracted waves retained. Square bullets: 
N e = N p evanescent diffracted waves retained. Results for the far-field are 
practically identical. 
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Figure 9: As in Fig. El with R = lOOAo, N e = N p , and various values of scale 
ratio v = difd\ and average spacing d av . (a): d av = 0.5Ao- (b): d av = 0.75Ao- 
Square bullets: v = 0.9. Circular bullets: v = 0.75. Triangular bullets: 
v = 1/r. 
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Figure 10: 101-element modified-Fibonacci array with d av = 0.874Ao and 
r] = 0. Radiation pattern (array factor) \F(k sm9)\ in ()10ajl (with \m\ < 50) 
for various values of scale ratio v, illustrating the controllable multibeam 

capability (see Sec. Q. v = 0.25 (S 01 = -1.83dB); v = 0.5 

(S" 01 = -6.37dB); v = 0.9 (S 01 = -23.3dB). Note that minor sidelobes 

never exceed the periodic-array sidelobe level (-13dB). 
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